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Abstract:  A method  for  diagnosing  multiple  element  defects  of  rolling  bearings  has  been  investi- 
gated. The  method  combines  the  time  synchronous  averaging  and  envelope  spectral  analysis 
techniques  to  produce  spectra  of  synchronously  averaged  envelope  signals  with  a range  of 
synchronous  frequencies.  The  spectra  are  displayed  in  the  domains  of  synchronous  period  versus 
frequency  to  result  in  the  sync-period  frequency  distribution.  The  distribution  separates  the  char- 
acteristic defect  frequencies  and  their  associated  sidebands  in  the  synchronous  period  axis.  This 
analysis  technique  makes  it  possible  to  detect  and  diagnose  multiple  defects  appearing  on  differ- 
ent elements  of  rolling  bearings.  Another  main  benefit  of  the  method  is  the  significant  noise 
reduction  by  both  enveloping  and  synchronous  averaging  process.  Analysis  results  from 
computer  simulated  data  and  experimental  simulated  data  are  presented. 
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1.  Introduction 

Time  synchronous  averaging  is  a powerful  tool  to  extract  periodic  waveforms  from  noisy  signals. 
It  has  been  used  extensively  in  diagnosing  gear  defects.  However,  this  technique  has  not  found 
wide  acceptance  for  diagnosing  defects  of  rolling  bearings,  because  of  the  asynchronous  rotation 
[ 1 ] between  the  shaft  and  the  rolling  elements  or  cage  of  the  rolling  bearing,  and  the  approxima- 
tion of  bearing  characteristic  defect  frequencies  [2],  On  the  other  hand,  the  high  frequency  reson- 
ance technique  (HFRT)  is  widely  accepted  in  bearing  defect  diagnosis.  For  single  defect  cases,  it 
has  proven  successful  under  reasonable  signal  to  noise  ratio  (SNR)  conditions,  provided  that  the 
shaft  is  running  at  a constant  speed  and  there  exists  no  skidding  of  the  rolling  elements.  But  this 
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technique  has  proven  less  than  satisfactory  in  detecting  bearing  defects  under  poor  SNR  condi- 
tions [3].  Also,  to  the  authors  knowledge,  the  successful  application  of  this  technique  to  the  diag- 
nosis of  multiple  element  defects  in  bearings  has  not  been  reported. 

In  diagnosing  multiple  faults,  Igarashi  and  Kato  [4]  investigated  the  vibration  of  a ball  bearing 
with  multiple  defects  on  the  race  surface  of  either  the  inner  or  outer  ring  under  given  thrust 
loads.  Using  an  FFT  analyzer  and  examining  the  effect  of  the  measuring  position  of  the  outer 
ring  on  the  vibration,  they  clarified  the  fundamental  characteristics  of  ball  bearings  with  multiple 
faults,  and  they  also  established  an  approach  to  locate  the  faults  and  indicate  the  number  and  size 
of  the  defects.  Mcfadden  and  Smith  [5]  described  an  investigation  of  a bearing  with  multiple 
point  defects  based  on  their  model  [6]  of  a single  point  defect.  The  vibration  envelope  spectra 
produced  by  a bearing  with  multiple  inner  race  defects  were  illustrated  by  the  reinforcement 
and  cancellation  of  spectral  lines  lue  to  the  different  phase  angle  of  the  vibrations  inaaced  by 
diffe:  ant  faults.  Braun  and  Datner  [7]  discussed  a bearing  fault  diagnostic  method  based  on  time 
averaging.  The  RMS  level  of  the  averaged  signal  was  utilized  to  indicate  the  relevant  defects 
according  to  their  characteristic  defect  frequencies.  This  method  effectively  extracts  the  periodic 
signals  associated  with  bearing  faults  from  noisy  structural  vibrations  and  is  less  sensitive  to 
structure  parameters  than  general  spectral  techniques.  However,  the  RMS  approach  seems  not 
capable  of  identifying  which  frequency  components  contribute  the  most  to  the  RMS  value, 
because  one  averaging  period  could  involve  several  frequency  components.  Moreover,  due  to  the 
low  energy  of  the  vibration  caused  by  bearing  faults,  the  characteristic  defect  frequencies  may  be 
very  difficult  to  be  discerned  when  the  shaft  rotation  related  frequency  components  (high  energy) 
are  closely  adjacent  to  those  defect  frequencies. 

The  high  frequency  resonance  technique  is  effective  in  excluding  the  shaft  rotation  related 
frequency  components,  but  it  has  no  advantage  in  separating  multiple  defect  frequencies  and 
eliminating  high  frequency  noise  disturbances.  In  cases  where  bearings  have  multiple  defects 
located  in  different  elements,  each  defect  may  have  different  shape  and  size.  Thus  each  defect 
may  excite  a different  resonance  (structural  and/or  transducer's)  to  a different  extent.  Using  the 
ordinary  high  frequency  resonance  technique,  only  one  resonance  is  isolated  by  a bandpass  filter 
for  envelope  detection.  Therefore,  where  multiple  faults  exist,  multiple  defpot  frequency  compo- 
nents are  likely  to  be  mixed  together,  provided  all  defects  excite  the  same  resonance  to  a compar- 
able extent.  Alternatively,  if  a particular  resonance  is  dominated  by  one  fault  excitation,  then 
only  that  fault  will  be  distinguishable. 

In  this  paper,  we  investigate  a technique  which  introduces  an  averaging  process  into  the  high 
frequency  resonance  technique  and  expands  the  bandwidth  of  the  bandpass  (or  even  high  pass) 
filter  to  include  several  high  frequency  resonances  rather  than  just  a single  isolated  resonance. 
The  time  synchronous  averaging  process  is  performed  after  the  envelope  detection.  Due  to  the 
uncertainty  of  the  bearing's  characteristic  defect  frequencies,  the  synchronous  frequency  is  swept 
around  the  nominal  defect  frequency  until  the  true  frequency  is  identified.  Spectral  analysis  is 
then  carried  out  for  the  resulting  averaged  envelope  signal.  The  spectra  are  finally  displayed  in 
the  so-called  syncperiod  - frequency  domain  to  form  the  syncperiod  - frequency  ( SPF)  distribu- 
tion. In  SPF  distributions,  the  individual  characteristic  defect  frequencies  and  their  associated 
modulation  sidebands  are  well  separated  along  the  syncperiod  axis.  The  exact  frequency  content 


can  then  be  acquired  so  as  to  diagnose  the  condition  of  the  bearing  with  multiple  element  defects. 
It  is  shown  that  the  enveloping  and  synchronous  averaging  processes  used  in  this  method  signifi- 
cantly improve  the  signal  to  noise  ratio. 

2.  Sync-Period  Frequency  (SPF)  Analysis 

In  this  section,  a method  which  combines  the  time  synchronous  averaging  technique  and  the  high 
frequency  resonance  technique  is  described.  From  the  high  frequency  resonance  technique,  we 
know  that  the  envelope  signal  generated  by  localized  bearing  defects  can  be  featured  as  a series 
of  exponentially  decayed  impulses  with  the  periods  related  to  various  defects  [6,  7],  For  the 
remaining  part  of  this  section,  we  are  only  dealing  with  the  envelope  signals. 

2.1.  Time  Synchronous  Averaging  of  a Bearing  Signal  and  Its  Spectrum 

Using  standard  envelope  detector,  the  raw  signal  is  firstly  bandpass  filtered  in  the  high  freouency 
resonance  region.  The  resulting  signal  consists  of  the  bearing  defect  induced  signal  embedded  in 
band-limited  Gaussian  noise.  Suppose  x(t)  is  such  an  envelope  signal,  and  it  contains  frequency 
components  corresponding  to  the  harmonics  of  the  characteristic  defect  frequencies  and  their 
associated  sidebands.  It  is  decomposed  into  A segments,  each  of  same  length  Ts.  Each  segment  is 
separated  by  a variable  period  X (X  > Ts  - ie.  no  overlap).  Summing  up  all  these  segments  yields: 

y(f,x)=  I x(t+nx),  te  [OX]  (1) 

n- 0 


where  t is  the  averaging  period,  Ts  is  the  segment  length.  According  to  the  time-shift  theorem  of 
Fourier  transform,  the  spectrum  of  (1)  Y(f,x)  will  be  [8,9] 
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the  amplitude  spectrum  will  then  yield 
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which  is  the  original  amplitude  spectrum  IX(/)I  multiplied  by  the  transfer  function  of  a comb  filter 
\H(f)\  = lsin(Arc/c)/sin(Jt/x)l.  By  examining  the  nature  of  the  function  \H(f)\  shown  in  Figure  1 [8, 
9],  we  can  see  that  a)  it  is  a periodic  function  with  the  period  of  /=  1/x;  b)  for  large  A the  func- 
tion appears  like  a train  of  sine  functions  and  it  becomes  an  impulse  train  when  A goes  to  infin- 
ity, with  period/=  1/x  ; c)  its  main  lobes  peak  amplitude  is  A,  they  are  located  at/=  n/x  ( n=0,l , 
...  ,N-1),  its  first  zero  is  at  / = 1/Ax  and  thus  the  main  lobes  have  the  width  of  B=2/Ax;  d)  its 
main  lobe  3 dB  point  is  at  0.45/Ax,  and  thus  3 dB  bandwidth  is  fi1dB  =0.9/Nx;  e)  its  side  lobe  peaks 
decrease  with  increasing  A,  the  maximum  side  lobe  peak  is  12^5  down  from  the  main  lobe  peak 
when  A=5,  and  13 dB  down  when  A=20,  it  approaches  the  limit  value  of  -\3.3dB  as  A increases. 
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Figure  1.  One  period  of  function  lsin(Mt/t)/sin(7t/t)l  [9] 

As  we  mentioned  early,  \X(f)\  has  the  frequency  components  of  mf„,  mf,  ± lrfr,  mfb  ± kfc  (m,k-0, 1 , 
2,...)  where  /,  /,  / denote  the  characteristic  defect  frequencies  of  outer  race,  inner  race  and 
ball/rollei  respectively,/  denotes  shaft  rotation  frequency,/  cage  frequency.  The  relative  ampli- 
tude of  these  components  depend  upon  the  damping  factor  of  the  bearing  system,  the  load 
distribution  and  the  variation  of  transmission  path  [1],  When  the  synchronous  period  t equals  an 
integer  multiple  M of  the  interval  corresponding  to  one  of  the  above  frequencies,  for  instance 
t -MTn-M/f0 , the  frequency  components  from  both  signal  and  noise  (except  / and  its  harmonics) 
are  rejected  and  a periodic  signature  with  a period  of  Ta  is  extracted  by  the  averaging  process  of 
(1)  (see  [10]).  Therefore,  the  spectrum  \Y(f,x)\  in  (3)  will  be  dominated  by/  and  its  harmonics. 

2.2,  Signal  Processing  Considerations 

Due  to  the  uncertainty  of  the  characteristic  defect  frequencies  of  the  bearings,  it  is  necessary  to 
identify  the  true  synchronous  period  by  searching  around  the  nominal  frequencies.  The  signal 
processing  scheme  that  is  used  to  find  the  true  synchronous  period  is  based  on  a search  strategy 
(Figure  2). 


Figure  2.  The  flowchart  of  Sync-Period  Frequency  Analysis 
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In  implementing  the  above  scheme  digitally,  several  factors  have  to  be  considered.  In  this 
section,  we  will  discuss  the  effect  of  sampling  rate/,  data  segment  length  Ts,  synchronous  period 
t,  the  number  of  averages  N and  the  search  spacing  At  on  the  performance  of  SPF  analysis.  For 
clarification,  we  tabulate  all  factors  concerned  in  SPF  analysis  as  following: 


fs 

sampling  frequency,  Atmn  = 1// 

T 

total  data  length 

Ts 

segmental  data  length 

B 

main  lobes  bandwidth  of  the  comb  filter,  B=2/Nx,  A/>  B 

fi3dB 

main  lobes  3dB  bandwidth  of  the  comb  filter,  Bmb  = 0.9/Nx 

f 

defect  frequency  components  to  be  detect 

A f 

frequency  accuracy,  defined  as  A/=  M[l/x-l/(x+Ax)],  (A/ = / Ax/x,  if  Ax  « x) 

df 

frequency  resolution  from  the  Fourier  transform,  df  = 7/T 

M 

integer  multiples  of  the  defect  periods 

X 

synchronous  period,  x = Mlf,  X ^ / (x  and  Ts  are  in  an  equivalent  order) 

Ax 

sync-period  search  spacing,  Ax/x  = Af/f  (ifAx«x),  Ax  jsl// 

N 

number  of  averages,  N = mod  (T,  x ) 

To  ensure  the  correct  detection  of  frequency  components,  we  are  particularly  concerning  about 
the  search  spacing  At  with  respect  to  the  frequency  accuracy  A/  (frequency  resolution  obtained 
from  the  syncperiod).  Large  Ax  may  fail  to  identify  the  bearing  defect  frequency,  whereas  too 
small  Ax  may  cause  a tremendous  increase  of  computing  time.  In  general,  Ax  is  determined  such 
that  the  frequency  accuracy  A/ is  larger  than  (or  equal  to)  the  bandwidth  [10]  of  comb  filter's 
main  lobes  shown  in  Figure  1.  According  to  the  requirement  of  the  SPF  analysis,  tactors  like/s, 
/,  Ax,  A/  aiid  B are  primarily  selected,  and  other  factors  are  determined  by  the  formula  in  the 
above  table.  Figure  3 depicts  how  non-dimensional  measures  A f/f,  Ax/x,  Ax/,  A/x,  A f/B,  N and  M 
are  related.  The  calculated  factors  could  be  verified  by  these  non-dimensional  graphs. 


Figure  3.  The  relationship  of  A f/f,  Ax/x,  Ax/,  A/c,  Af/B,  N and  M.  (a)  A f/f  vs.  Ax/(or  xA f),  both  in  logarithm  scale, 
Af/f  = Ax/ /M-  (b)  Af/B  vs.  Ax/  (or  xA/),  both  in  logarithm  scale,  Af/B  = AxfN/2. 

For  example,  assuming  the  sampling  rate  is  100kHz  which  leads  to  a Ax^  of  0.01ms,  if  we  select 
Ax  = 0.04ms  and  A f—  0.25Hz,  then  to  detect  a frequency  component  of  250Flz  would  require  M 
to  be  10,  and  x to  be  40ms.  Furthermore,  if  we  select  B = A/=  0.25Hz,  N would  be  200  and  T 
would  be  8 s (800000  samples). 
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In  detecting  multiple  frequency  components,  we  need  to  sweep  the  synchronous  period  x in  a 
wide  range.  We  can  keep  Ax  unchanged  during  the  sweeping  while  using  different  M 's  to  associ- 
ate with  multiple  / 's.  In  the  above  example,  if  there  are  three  frequency  components  to  be 
detected,  eg.  152Hz,  198Hz,  247Hz,  M's  would  be  4,  6 and  10  respectively,  and  the  sweeping 
range  of  x would  be  26.32ms  ~ 40.5ms.  To  ensure  a successful  search,  the  range  could  be 
expanded  to  24ms  ~ 43ms.  If  a fixed  averaging  number  N (eg.  200)  is  selected,  the  frequency 
bandwidth  B will  vary  from  0.42Hz  to  0.23Hz,  and  Af  from  0.23Hz  to  0.24Hz.  Provided  the 
range  of  Af/B  is  acceptable,  fixed  N is  more  convenient  for  sweeping.  In  this  regard,  the 
frequency  resolution  for  low  frequency  components  will  be  poorer  due  to  the  smaller  Af/B. 

Following  the  averaging  process,  the  frequency  spectrum  of  the  averaged  envelope  signal  is  esti- 
mated. The  FFT  associated  frequency  resolution  of  the  spectrum  df  i*  therein  another  concern. 
However,  during  sweep  synchronous  averaging  the  sync-period  x ( =M/f)  can  give  a much  higher 
frequency  resolution  than  what  could  be  achieved  through  spectral  analysis.  The  required  df 
should  be  high  enough  to  allow  individual  frequency  components  to  be  distinguished.  For  the 
above  example,  the  spectral  resolution  df  is  about  38Hz  ~ 25Hz  (1/x). 

3.  SPF  Pattern  of  Multiple  Bearing  Defects 

The  SPF  patterns  provide  a graphical  exhibition  of  the  frequency  composition  of  a signal  in  the 
frequency-syncperiod  domain.  In  this  section,  we  present  a SPF  pattern  for  a combination  of 
different  types  of  element  defects  using  nominal  characteristic  defects  frequencies  [7,  11]  of  a 
bearing.  The  main  uncertainty  in  estimating  the  characteristic  frequencies  is  the  contact  angle 
between  rolling  elements  of  the  races,  which  is  dependent  on  the  type  of  bearing  and  loading 
conditions.  Another  source  of  uncertainty  is  the  diameter  variations  of  the  balls. 

In  the  cases  where  multiple  defects  in  different  bearing  elements  exist  and/or  the  SNR  is  poor, 
envelope  SDectra  [6]  prove  unsatisfactory  for  the  identification  of  individual  bearing  defect 
frequency  omponents.  SPF  pattc  ns  indicate  where  the  characteristic  defect  frequencies  and 
their  sidebands  are  located  in  a two  dimensional  domain  (notice  that  the  relative  intensity  is  not 
shown),  referred  here  as  'syncperoid  - frequency  domain'.  Syncperiod  axis  represents  the  sweep- 
ing range  for  the  synchronous  period  x which  is  an  integer  multiple  of  the  characteristic  defect 
periods,  the  frequency  axis  is  used  to  define  the  frequency  contents  corresponding  to  a given  the 
syncperiod.  In  fact,  with  small  Ax  the  syncperiod  x contains  much  finer  frequency  resolution  than 
frequency  axis.  However,  without  the  frequency  axis,  the  syncperiod  along  would  not  be  able  to 
distinguish  which  frequency  components  this  period  multiple  ( x=M/f ) is  associated  with.  For 
instance,  a syncperiod  of  24ms  could  be  either  6 times  4ms  (250Hz)  or  4 times  6ms  (166.67Hz). 

Now  we  consider  a rolling  bearing  (NSK  EN202)  with  different  types  of  defects  and  present  the 
SPF  pattern  associated  with  a combination  of  various  defects.  This  type  of  bearing  is  deliberately 
chosen  for  the  near  harmonics  of  various  defect  frequencies,  ie.  3 fr  and /,/-/  and/,  2/  and  f +/ 
etc.  are  closely  adjacent,  which  makes  the  diagnosis  of  defects  more  complex.  Dimensions  and 
defect  frequencies  of  the  bearing  are  listed  below.  The  SPF  pattern  is  shown  in  Figure  4. 
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Figure  4.  SPF  pattern  for  a combination  of  outer  race,  inner  race  and  roller/ball  defects  ('o'  - outer  race  defect  related 
components;  - inner  race  defect  related  components; '+'  - ball/roller  defect  related  components;  all  arrows  are 

pointing  at  the  fundamental  defect  frequencies) 

4.  Diagnosis  of  Multiple  Defects  Located  in  Different  Bearing  Elements  Using  SPF  Analysis 

In  this  section,  we  will  discuss  the  application  of  the  SPF  analysis  to  detect  multiple  defects  in 
different  bearing  elements.  Results  from  both  computer  synthesised  and  experimental  simulated 
data  are  described. 

4.1.  Diagnosis  of  Multiple  Defects  Generated  by  Models 

The  synthetic  vibration  signal  used  in  this  section  is  a bearing  signal  generated  via  Braun’s  model 
[2,  7]  plus  other  vibration  sources  (eg.  unbalanced  shaft,  etc.)  and  noise.  The  following  example 
presents  the  SPF  distribution  from  synthetic  signal  generated  by  the  model.  It  consists  of  an 
impulse  series  induced  by  all  three  types  of  defects,  shaft  rotation  related  vibrations  (narrow  band 
but  high  power)  and  broadband  white  noise.  The  maximum  impact  amplitudes  are  35,  90,  100 
for  the  outer  race  defect,  inner  race  defect  and  ball  defects  respectively.  The  noise  standard  devi- 
ation <yn  is  30.  It  is  also  assumed  r our  synthetic  data  that  three  different  types  of  defects  excite 
three  different  ..sonants,  ie.  /,=/5070Hz,  /2=1 8003Hz,  f3=2 11 10Hz,  and  with  different  time 
constants  of  exponential  decay,  a,=l/15,  02=03=2/15,  respectively.  Using  the  SPF  analysis 
described  in  section  2,  the  SPF  distribution  of  the  synthetic  signal  plus  broadband  white  noise  is 
generated  and  displayed  in  Figure  5.  All  factors  used  are  tabulated  below: 
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Comparing  this  SPF  distribution  with  the  SPF  pattern  shown  in  Figure  4,  it  is  found  that  the 
characteristic  defect  frequency  components  and  their  associated  modulation  sidebands  are  well 
separated  in  syneperiod  - frequency  domain.  For  a given  synchronous  period,  only  one  frequency 
component  and  its  harmonics  are  exposed.  It  can  be  seen  from  the  SPF  distribution  that  the  lower 
frequency  components  are  more  widely  spreaded  over  the  syneperiod  domain.  The  noise  level 


has  been  significantly  attenuated  through  enveloping  (3dB  improvement)  [12]  and  averaging  (a 
factor  of  Jn  ) [10,13], 

Sync. Period-Frequency  Distribution  (25dB  down) 


Figure  5.  SPF  distributions  from  synthetic  signal  plus  broadband  noise  with  O„=30 

In  fact,  the  original  signal  for  Figure  5 has  a very  poor  signal  to  noise  ratio  (-15.5  dB).  From  time 
domain,  there  is  no  evidence  of  any  characteristic  defect  frequency  related  components  because 
of  high  level  noise  disturbance  and  the  intertwinement  of  multiple  frequency  components.  Using 
the  envelope  spectral  technique  (Figure  6),  it  is  also  very  difficult  to  identify  individual  char- 
acteristic defect  frequencies. 


Figure  6.  (a)  Original  time  signal  for  Figure  5 and  (b)  the  envelope  spectrum  (1024  averages)  of  (a) 

4.2.  Diagnosis  of  Multiple  Defects  With  Experimental  Data 

Two  NSK  EN202  bearings  were  used  for  this  experimental  investigation.  The  artificial  defects 
were  implanted  by  either  nitric  acid  corrosion  or  electric  discharge  erode.  Bearing  #1  has  one 
defect  on  both  outer  race  and  inner  race  (nitric  acid  corrosion  - about  2-3  mm  in  diameter, 
depth  unknown  but  very  shallow).  Bearing  #2  has  one  ball  defect  implanted  (electric  erode  - 
diameter  0.38mm,  depth  1mm),  and  its  cage  was  dismantled  and  re-riveted  after  the  introduction 
of  the  defect.  The  analysis  results  are  shown  in  Figure  7 and  8.  All  SPF  analysis  parameters  are 
identical  to  the  examples  in  section  4.1  except  the  range  of  syncperiod  x for  bearing  #2  (from 
99.01-107  ms).  The  experimental  conditions  are  listed  in  the  following  table: 


shaft  speed 


transducer 


B&K  8307  (#1) 
B&K  4393  (#2) 


load  applied  sampling  rate  /„,  /,,  fc 


240N  radial  (#1)  100  kHz  152.47,  247.53,  198.52,  19.06  Hz 

260N  radial  (#2)  (nominal  defect  frequencies) 


Figure  7. 
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Figure  8.  SPF  distribution  of  Bearing  #2,  the  signal  was  bandpassed  at  40~43kHz  before  enveloping 

5.  Discussions  and  conclusion 

By  baseband  spectrum  analysis  of  vibration  signal  of  bearing  #1,  it  is  found  that  a series  of  high 
frequency  resonances  have  been  excited.  Envelope  demodulated  spectra  at  several  resonances 
show  both  inner  race  and  outer  race  defective  frequencies  and  associated  shaft  rotation 
frequency,  but  all  those  frequency  components  are  intertwined  together.  The  SPF  distribution 
shown  in  Figure  7 separates  those  frequencies  along  the  syncperiod  axis  and  shows  the  defective 
frequencies  and  the  modulation  sidebands  explicitly.  The  outer  race  defect  frequency  is  derived 
from  the  point  at  105.23ms  on  syncperiod  axis  and  about  150Hz  on  frequency  axis,  which  corre- 
sponds to  152.05Hz  (M  is  16  here).  The  inner  race  defect  frequency  (248.43Hz,  M is  25)  is 
corresponding  to  the  point  at  100.63ms  and  250Hz.  The  SPF  distribution  of  bearing  #2  displays  a 
multiple  defect  pattern  rather  than  a single  ball  defect  pattern.  The  indication  of  inner  race  and 


outer  race  defects  are  very  strong.  This  curious  phenomenon  suggests  that  inner  race  and  outer 
race  defects  may  have  been  produced  during  dismantling  and  re-assembly  of  the  bearing.  The 
envelope  spectrum  from  a low  frequency  resonance  (2kHz  ~ 4kHz)  reveals  inner  race  fault 
related  components.  The  outer  race  defect  frequency  dominates  the  envelope  spectrum  demodu- 
lated about  a 26k~28kHz  resonance.  Ball  defect  components  could  be  seen  from  the  envelope 
spectra  of  2kHz  highpassed  and  17k~23kHz  bandpassed  signals.  The  SPF  distribution  from  the 
40~43kHz  bandpassed  signal  seems  to  give  the  best  presentation  of  all  three  types  of  defects.  By 
inspecting  Figure  8 closely,  the  ball  defect  frequency  and  its  sidebands  are  only  discernible  with 
ease  around  100msec  region  due  to  the  small  size  of  the  defect.  The  cage  component  at  around 
104.65ms,  which  corresponds  to  19.1 1 Hz  (M=2),  is  an  important  indication  of  ball  damage. 

In  Figure  5,7,8,  we  notice  that  the  low  frequency  components  tend  to  spread  over  a wide  syncper- 
iod  region.  The  reason  for  that  may  be  explained  by  uncertainty  principle  [13].  Because  the  sync- 
period  domain  is  virtually  a frequency  scale,  a certain  length  of  data  segment  Contains  less  cycles 
for  low  frequencies  than  those  for  high  frequencies,  which  results  in  poorer  frequency  resolution. 

The  above  discussion  suggests  that  the  syncperiod  frequency  analysis  combines  the  benefits  of 
both  high  frequency  resonance  technique  and  time  synchronous  averaging  technique  in  overcom- 
ing problems  like  low-duty  bearing  related  vibration  and  asynchronous  rotation  of  rolling 
element  bearings.  It  is  therefore  concluded  that  the  syncperiod  frequency  analysis  is  potentially 
an  effective  method  for  the  diagnostics  of  multiple  defects  presented  in  different  bearing 
elements. 
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